Non-local order in gapless systems: Entanglement Spectrum in Spin Chains 
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We show that the entanglement spectrum can be used to define non-local order in gapless spin 
systems. We find a gap that fully separates a series of generic, high 'entanglement energy' levels, 
from a flat band of levels with specific multiplicities that uniquely define the ground-state, and 
remains finite in the thermodynamic limit. We pick the appropriate set of quantum numbers, and 
then partition the system in this space. This partition corresponds to a very non-local real-space cut. 
Despite the fact that the Laughlin state is bulk gapped while the antiferromagnetic spin chain state 
is bulk gapless, we show that the S = i Heisenberg antiferromagnet in one dimension has an en- 
tanglement spectrum almost identical to that of the Laughlin Fractional Quantum Hall state in two 
dimensions, revealing the similar field theory of their low-energy edge and bulk excitations respec- 
tively. In addition, we investigate the dimerization transition from the perspective of entanglement 
gap scaling. 

PACS numbers: 73.43.2013f, 11.25.Hf 
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Gapped topological phases of matter usually lack local 
order parameters that can distinguish them from trivial 
ones. In the presence of a gap in the excitation spectrum, 
several non-local indicators of the topological nature of 
the topological phase, such as ground-state degeneracy 
on compact high-genus manifolds [H as well as the struc- 
ture of edge modes and their scaling exponents Q exist, 
but do not fully describe the topological phase. In the 
absence of an excitation gap, the nature of quantum or- 
der is much less understood, especially in systems lacking 
a local order parameter. For example, torus degeneracy 
cannot distinguish gapless phases, since the mixing of 
infinitesimally low-energy excitations can lead to many 
new states quasi-degenerate with the true ground-state. 
There are several examples of gapless phases of matter, 
most notably in one-dimension (but also in higher dimen- 
sions [3]), which are known to display highly non-trivial 
properties. The parade example is the spin-^ Heisen- 
berg chain, a spin liquid which exhibits new elementary 
excitations, the spinons [H, HI, which fractionalize out of 
the usual antiferromagnetic spin wave. While for gapped 
systems such as the AKLT spin chain, non-local order 
parameters have been found, no similar quantity exists 
for gapless spin-disordered chains. 

The advent of quantum information theory has led 
to new concepts, such as quantum entanglement, which 
have proved useful in further characterizing quantum or- 
der. The majority of research to date has concentrated 
on entanglement entropy of a state \ip ). Specifically, let 
the Hilbert space for an A-body system be written as 
a direct product H = H A ® rl B . It is easiest to imag- 
ine real space partitions, where the division corresponds 
to some prescribed real space 'cut'. Then p A = Tr B (p), 
where p — \ipo)(4>o\i ^ s the reduced density matrix for 
the A component. The quantity S A — —Tr A (p A hip A ) 
is the entanglement entropy, which provides a measure 
of the quantum entanglement of |^ Q ) with respect to the 



partitioning (A,B). 

Different scalings of this quantity with the size of A 
for gapped and gapless systems have been proposed and 
proved @. Still, the entanglement entropy remains a 
single number, and can only provide a limited character- 
ization of topological order. To this end, Li and Haldane 
(LH) in a recent insightful paper [7] proposed and nu- 
merically substantiated that the entanglement spectrum, 
i.e. the full set of eigenvalues of p Al partitioned with 
respect to the equator of the quantum Hall sphere 
provides a near complete picture of the topological order 
in v = I fractional quantum Hall effect (FQHE) states. 
Writing the eigenvalues of p A as e~^< , where ^ is an en- 
tanglement level, LH showed that the low level spectrum 
for generic gapped v = | states exhibited a universal 
structure, related to conformal field theory, and sepa- 
rated from a non-universal high energy spectrum by an 
entanglement gap which was finite in the thermodynamic 
limit. This gap itself was proposed as a 'fingerprint' of 
the topological order present. It was subsequently shown 
that the entanglement spectrum can meaningfully distin- 
guish among states which have similar finite size overlap 
with |^ ) but different edge structures [9(. Other bulk- 
gapped phases such as topological insulators have also 
been shown to reflect their nontrivial nature in their en- 
tanglement spectra HI II |- 



In this Letter, we show how to use the entanglement 
spectrum to characterize non-local order in gapless sys- 
tems and phases of matter, with emphasis on spin-i spin 
chains. We find that simple real space partitions fail 
to reveal the nontrivial fundamental structure in these 
phases. Rather, we find the cut must be made in Fourier 
(momentum) space. The local cut in momentum space 
corresponds to a highly non-local cut in real space. We 
only use the ground-state wavefunction of the system and 
do not assume any knowledge regarding the gapless ex- 
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citation spectrum. If the system is in the SU (2) x Wess 
Zumino Witten (WZW) universality class for spin- 5, we 
show that the low-entanglement level part of the entan- 
glement spectrum has a special form similar to that of 
the Laughlin [12J FQHE state. The high entanglement 
levels are separated by a large entanglement gap from the 
low entanglement WZW portion of the spectrum. This 
recapitulates the situation for the ground state of the 
Coulomb Hamiltonian on the quantum Hall sphere Q, 
but with two important differences: (i) we are here deal- 
ing with a spin system, and (ii) we find a complete entan- 
glement gap at all values of the total momentum, whereas 
LH found a gap only for a finite range of angular momen- 
tum L z . We numerically show that the entanglement 
gap is finite in the thermodynamic limit for the Heisen- 
berg chain. As one result of our work, one can classify 
whether a generic ground-state spin- \ wavefunction is in 
the WZW universality class, without resorting to often 
unreliable computations of critical exponents. We also 
show that the entanglement spectrum is sensitive withe 
respect to the dimerization transition - the flow of entan- 
glement levels matches the field theory prediction for the 
canceling of a marginal operator at exactly the dimeriza- 
tion transition. 

We consider spin-i spin chains on even-membered 
rings with periodic boundary conditions (PBCs). The N 
sites are placed on a circle of radius unity and are hence 
described by the iV th roots of unity: Zj — e 2mj / N ; j £ 
{1, . . . , N} (see Fig. [TJ. Without loss of generality, any 
ground state wavefunction of an SU(2)-invariant Hamil- 
tonian can be written in the form 
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where F ) = 1 1 ' ' 4 t ) is the ferromagnetic state. The 
sum extends over all ways to distribute the K = 
down-spins on the ring, and the weights ip( z j 1 y ■ ■ ■ > z j ) 
depend only on the position of the spin 4. sites (Holstein- 
Primakoff representation). We further presume that |^ ) 
is a translationally invariant total spin singlet (otherwise, 
in most cases, the state could already be characterized by 
some other local order parameter). In particular, these 
are also the conditions for spin liquid states. Mathemat- 
ically stated, we require J2j Sj~\^o) = an d 
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where Q = or 7r is the total crystal momentum, depend- 
ing on whether K is even or odd, respectively. Note that 
the last condition is a generic property independent of 
the particular ground-state wavefunction studied. Spin 
rotational invariance also implies that the wavefunction 
is invariant under interchange of t and 4- coordinates. 

Typically, calculations of entanglement and entangle- 
ment entropy proceed by cutting the system spatially into 



Zi Zj 






I>\ 2 




\ 1 




IN 



14 _ 

12 _ 
10 _ 
8 _ 

6 
4 

2- 



1 



Sa 



FIG. 1: (color online) Left: Unit circle description of a spin 
chain with PBCs; the sites lie at the N th roots of unity. The 
chord distance specifies the spatial separation of any two sites 
Zi and Zj. Right: Typical low- level entanglement spectrum for 
a spatially bisected N = 30 Heisenberg chain (ground state). 
Entanglement levels along the vertical axis are plotted versus 
total spin Sa quantum number. Each level has a 2Sa + 1 
multiplet degeneracy. No immediate structure is visible. 



two disconnected parts and then computing some com- 
bination of the eigenvalues of the reduced density matrix 
which decomposes in independent sectors indexed by the 
values of the total S z . If we choose to perform this cut, 
compute the eigenvalues of the density matrix, and plot 
them versus the good multiplet quantum number Sa, 
one observes an unremarkable structure such as shown 
in Fig. [T] To obtain relevant information, one must then 
analyze the distribution of eigenvalues of the entangle- 
ment spectrum [13| or perform calculations of the entan- 
glement entropy [g, |l4 1 . This suggests that a spatial cut 
is not the most useful way to understand the underly- 
ing structure of the spin chain, as no differences in the 
systems can be suitably resolved [131 ]. 

However, real space coordinates are only one option of 
a set of good quantum numbers in the system. We can 
also Fourier transform the spin operator on each site and 
go to a magnon-type description. We define the Fourier 



transform of the spin flip operators S m = jfYl 
The ground state is then written as 
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|* }= E 4>(rn 1 ,...,m K )S- i ---S- K |F> , (3) 



m 1 ,...,m K 



where 



^(m 1; . . .,m K ) = N 



- N - K E *r 



(4) 

We stress that while the momenta to, also belong to the 
set m 6 {1, . . . , N}, the momentum space spin flips do 
not behave as hard core bosons. 

The non-orthogonal momentum state basis can thus 
be represented by bosonic occupation numbers n m for 
crystal momentum m. The total particle number is 
K = E m n m, and e'Q = H m e 27rmn ™ /N = e 2mM / N ', 
where < n m < K V m. For generic wavefunctions, 
there exist multiple total momentum sectors separated by 
27r, although special spin-chain wavefunctions can have 
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FIG. 2: (color online) a) Entanglement spectrum of the Heisenberg model for N=22 sites, with a cut region containing 5 
magnons. The eigenvalues £ are plotted versus the total momentum subspaces of region A, Ma- One observes a low level 
part well separated from higher level contributions. This entanglement gap persists in the thermodynamic limit as shown in 
the upper inset. For the Haldane Shastry spin chain, exactly the same low-lying levels are found, where the higher levels are 
shifted to infinity, b) Entanglement spectrum of the K = 11 particle v = 1/2 FQH Coulomb state for a half cut on the sphere 
containing 5 particles. The counting of the low level states is exactly the same. The lower insets in (a) and (b) show the 
analogous level counting 1, 1, 2, 3, 5, 7 for the two states. 



compact support, i.e. weight in only one of these sectors. 
We normalize the states by the bosonic normalization of 
the magnon basis. Now consider a cut in momentum 
along the middle of the momentum orbitals. Both sepa- 
rate regions A and B can be decomposed with respect to 
number of particles and total momentum, which is sub- 
ject to the constraint Na + Nb = K and Ma + Mb = M. 

Let us first consider the ground state of the S = h 
Heisenberg model, H = ^ . Sj ■ Sj+i, for N = 22 sites. 
For systems of this size, the explicit Fourier transforma- 
tion Eq. H] exceeds numerical limits. Instead, we diago- 
nalize the spin Hamiltonian directly in momentum space, 
with matrices containing up to several billions of scatter- 
ing elements. For the Heisenberg ground state, we find 
that 98% of the weight of the state is present in the sector 
with M — K 2 . In the thermodynamic limit, we conjec- 
ture that the same low-energy structure of the entangle- 
ment spectrum will be present in all momentum sectors. 
The decomposition of the Heisenberg ground state in mo- 
mentum notation for a middle cut is presented in Fig. [3] 
We first observe a low entanglement level part clearly 
separated from levels higher up. For the low-lying states, 
from right to left, we observe the counting 1,1,2,3,5,7 
of levels of a U(l) boson mode until we hit a limit due 
to finite size of the system. This counting is equivalent 
to that for the excitations of the low energy field theory 
for the spin-i Heisenberg chain (the SU(2)i field theory is 
equivalent to that of a free U(l) boson, whose level count- 
ing is equal to the number of partitions of its momentum 
above the ground-state) . We then compare to the entan- 
glement spectrum of the Coulomb FQHE state on the 
sphere with the same number of particles. The L A quan- 
tum number on the quantum Hall sphere corresponds to 
the total momentum Ma for the cut region, where the 



convention is that the momentum counting starts from 
the cut, i.e. for maximum Ma, the occupied momentum 
orbitals have maximal orbital distance to the cut. This 
way, equal number of particles beyond the cut and equal 
absolute values of Ma and L A specify the same configu- 
rations in bosonic notation. For the Heisenberg ground 
state on the chain and the Coulomb state on the sphere, 
the entanglement spectra are similar both by level se- 
quence as well as by state counting, except that in the 
spin chain case, a full entanglement gap can be defined 
between all the SU(2) 1 low entanglement levels and the 
higher lying ones (Fig. [5]). A finite size scaling analysis 
of the entanglement gap for the Heisenberg ground state 
is shown in the inset of Fig. [2^.. Despite the relatively 
small number of sites that are available to us on small 
computers, we are confident that the entanglement gap 
remains finite in the thermodynamic limit. 

The entanglement property appears to be intimately 
linked to the low-energy properties of the system. Other 
regimes of the spin-i system possess different signa- 
tures of their entanglement spectrum. Adding a second- 
neighbor interaction, H = J2j Sj ■ Sj+i + A Sj ■ Sj + 2, we 
find that as we shift from the Heisenberg point (A = 0) 
to the Majumdar-Ghosh [l5[ point (A = |), where the 
ground state is a dimer crystal, the entanglement spec- 
trum significantly rearranges, reflecting the failure of the 
gapless U(l) boson description of the problem in the 
dimerized phase. From the Sine- Gordon description of 
the model [l6[ that we know that the leading logarith- 
mic operator corrections vanish at the dimerization point. 
This corresponds to the vanishing of the Sine-Gordon co- 
sine term at A c . At this critical point, only subleading 
operators remain, and the theory should look more like 
that of a free-boson than even the Heisenberg point. We 
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find this is manifested in the entanglement spectrum. 
Upon increasing frustration from the Heisenberg point 
A = 0, we find that the lowest (i.e. leading correction) 
high-energy levels first move upwards, which yields an in- 
crease of the entanglement gap A (Fig. [3^). The dimcr- 
ization transition is given by the vanishing of the leading 
marginal operator and hence by the maximum entan- 
glement gap. Under this hypothesis, we obtain a value 
fo J2/J1 = 0.2f5 in the thermodynamic limit. This is 
close to the established 0.24 value and even more note- 
worthy as our largest size system is 22 sites. As the 
coupling of the Sine-Gordon cosine term changes sign at 
A c , the term is responsible to generate the dimerization 
gap. This is seen in the entanglement spectrum by the 
reduction of the entanglement gap, as the dimerization 
constant A is tuned across its critical value. We remark 
that our numerical effort was heavy: brute-force Fourricr 
transform of the wavefunction was possible only up to 18 
sites. We then re-formulated our hamiltonian in Holstein- 
Primakoff basis, linearized, and introduced a term that 
gives large energy for double boson occupancy on site. 
We were then able to reach 22 sites, with non-sparce ma- 
trices of several billion matrix elements. 

Beyond the same state counting of the universal low 
level portions of their entanglement spectra, there are 
several other analogies between the Heisenberg and 
Laughlin states. As was first discovered for quantum 
Hall trial states 17[, we likewise find for the Heisenberg 
state that, to a large extent, the weights of the differ- 
ent basis states obey a "product rule" in the sense that 
the weights of a Heisenberg state of size N can be com- 
posed out of weight products of sub-parts of the bosonic 
momentum state configuration {m 11 . . .m K }, given by 
Heisenberg states of accordingly smaller system size. 

This strong structural similarity between the Heisen- 
berg state and Laughlin state is not coincidence. The 
Heisenberg chain is in the same universality class as the 
Haldane-Shastry (HS) chain [HI, EH, whose ground state 
wavefunction is *(z 1; ...,z K ) = W^Azi - Zj) 2 n t Zi- 
Though of no consequence in real space (where the entan- 
glement spectra for HS are unremarkable) , we found that 
the lattice Fourier transformation over all sites (Eq. |4]) of 
the HS state yields exactly the weights of the monomials 
m '{n K ,n K in the if -particle Laughlin state. Hence 

the entanglement spectrum of the HS ground state, using 
a momentum cut, is identical in its counting to the en- 
tanglement spectrum of the Laughlin state, as shown in 
Fig. [2] 21]. We stress that while the entanglement levels 
will be different due to the different normalizations for 
the spin chain and the QH sphere, the counting is iden- 
tical. The HS state is known to be the conformal gaplcss 
phase of all spin models within the the SU (2)2^ WZW 
class, with special quantum group symmetries. In terms 
of the fractionalized spinon excitations, it is the model 
in which the spinons are free in the sense that they in- 
teract only through their fractional statistics. Toward 
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FIG. 3: (color online) a) Entanglement gap A versus dimer- 
ization coupling A for different system sizes. As we move 
away from the Heisenberg point A = 0, the entanglement gap 
rise, till it finally shrinks again within the dimerized phase, 
b) The dimerization transition is estimated from finite size by 
the 1/N scaling of the A with maximum entanglement gap for 
different system sizes. The extrapolation gives A£° ~ 0.215, 
i.e. in range of 10% of the exact result from CFT. 



the Heisenberg point, the spinons get dressed, but retain 
their fractionality. The entanglement spectrum of the HS 
state consists purely of the universal low spectral levels 
of the Heisenberg phase (all the other levels are shifted 
to infinity in £, corresponding to a largely degenerate O 
eigenvalue level of pa)- From there, the entering loga- 
rithmic CFT corrections appearing as finite high-energy 
entanglement levels can be nicely observed upon inter- 
polating between the ground states of the HS and the 
Heisenberg Hamiltonian. This underscores the fact that 
the HS and Heisenberg model belong to the same univer- 
sality class, with the same relation between the Laughlin 
state and the Coulomb state in the FQHE. 

We have introduced a suitable basis to fruitfully discuss 
entanglement spectra of quantum spin chains. Specifi- 
cally, we have analyzed the S = h Heisenberg ground 
state wavefunction and observe a universal set of states 
separated from higher levels by an entanglement gap per- 
sisting in the thermodynamic limit, bearing strong simi- 
larities to the Laughlin and Coulomb states in the FQHE. 
It is hence possible to determine the non-local structure 
of a spin system just by examining the ground state, 
without any reference to the concrete Hamiltonian, as 
we demonstrated for the Heisenberg state in its relation 
to the Haldane-Shastry state, as well the dimerization 
transition from the Heisenberg to the Majumdar-Ghosh 
state. From the present result and the FQH studies, a 
unified picture emerges: the low level structure of the 
entanglement spectrum reflects the properties of the low- 
energy excitations of the system; it is in a one-to-one cor- 
respondence with the elementary excitations of the cor- 
responding field theory. For FQH states, which are bulk 
gapped, an orbital cut as performed in 0, Q is similar 
to a spatial cut due to the localized nature of the Lan- 
dau orbitals on the sphere. The entanglement spectrum 
then reveals the nature of the low-energy excitations of 
the gapless edge system. For Heisenberg spin chains, the 
bulk is gapless, but the momentum cut does not longer 
correspond to a spatial cut, and the entanglement spec- 



trum reveals the nature of the bulk gapless excitations. 
In both cases, these low-energy excitations are described 
by U(l) bosonic field theory. 
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FIG. 4: Entanglement gap evolution of the interpolation 
Hamiltonian ([5]) with TV = 16 sites for the values 7 = 0.75 
and 1.50. Note that high entanglement levels start to appear 
immediately away from the HS point, resulting in an entan- 
glement gap. Already at 7 = 1.5, the spectrum already closely 
resembles the Heisenberg spectrum. 
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FIG. 5: Entanglement spectra of the v = i Laughlin (left) 
and HS (right) states for systems with 9 particles (or down- 
spins). The spectra consist solely of low level states; the en- 
tanglement gap is infinity. The edge mode counting exactly 
matches the prediction from the CFT. Both spectra possess 
exactly the same structure, while the HS spectrum - in abso- 
lute values of entanglement levels - is much flatter. 



Supplementary material: 

(i) Logarithmic corrections from the Haldane-Shastry 
point. We show how the entanglement spectrum changes 
from the Haldane-Shastry point to the Heisenberg point. 
This is best accomplished by the interpolation Hamilto- 
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l<3 



2simi (7|z, - Zj\/2) 



(5) 



where the Haldane-Shastry limit is reached for 7 — ► 
0. From the viewpoint of CFT, as discussed in the 
manuscript, the change in the entanglement spectrum 
relates to the logarithmic corrections entering from the 
conformally critical HS point. We find that the states 
enter in a flat band structured fashion corresponding to 
different orders of corrections, as shown in Fig |H 

(ii) Correspondence of Haldane-Shastry and Laughlin 



spectrum. As discussed in the manuscript, the Haldane- 
Shastry state and the Laughlin state basically have the 
same entanglement spectrum, which only differs in terms 
of the manifold where the states are defined. Accord- 
ingly, the state counting exactly matches. However, the 
Laughlin levels spread over a much bigger range of entan- 
glement energy, while the Haldane Shastry levels form an 
extremely flat band (Fig. [5]). 

(Hi) Convergence of Lanczos algorithm. The biggest 
system size are computed by diagonalizing the spin 
Hamiltonian in a hardcore boson model in momentum 
space. The distribution of off-diagonal elements imposes 
a considerable challenge for Lanczos optimization con- 
vergence to generate the eigenstate. Exploiting all sym- 
metries available in momentum space, the biggest sizes 
considered (N — 22) are matrices of dimension 3848879 
with O(10 9 ) scattering elements. Convergence within a 
standard error of 0(1O~ 8 ) for the Lanczos regeneration 
of the eigenstate demanded up to 50000 iterations. 

(iv) Fitting procedures for the entanglement gap. There 
are slightly different ways to define the entanglement gap, 
all of which yield very similar results. The convention 
used in the paper is 



A := Min„ 



.E{M< 



]{Min 5 {C ncric }-Max 5 {e 



universal 



}}• 



(6) 

i.e. the minimal gap between the highest low-energy uni- 
versal level and lowest high-energy generic level of all 
different momentum sectors of M A . Alternatively, one 
can define the entanglement gap globally over all different 
momentum sectors Ai := Min{£s cnoric }-Max{£ univorsa1 }. 
The latter Cciii be further generalized to ct A n where 
instead a just taking the lowest generic level, we average 
over the lowest n generic levels. If there were a spectral 
flow of only a single mode closes the gap, the difference 
in Ai and A n would be able to detect it. For the cases 
considered in the manuscript, however, we only observe 
a collective rise or lowering of levels sets extending over 
all M A subsectors, which is why all these gap definition 
nearly conincide in their result. 

(v) Spectral flow for the dimerization transition. A 
movie (movie-jlj2.gif) of the spectral flow from the 
Heisenberg point to the Majumdar Gosh point is added 
to the arxiv source files. 



This figure "movie-jlj2.gif" is available in "gif" format from: 



http://arxiv.org/ps/0912.0028v2 



